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ABSTRACT 


The power transfer efficiency between two circular apertures at an arbitrary 
range is obtained numerically. Two apertures can have generally different sizes 
and arbitrary taper illuminations. The effect of distance and taper illumination 
on the transmission efficiency are investigated for equal size apertures. The 
result shows that microwave beam power is more effective at close ranges, namely 
distances less than 2D 2 /lambda. Also shown was the power transfer efficiency 
increase with taper illumination for close range distances. A computer program 
was developed for calculating the power transfer efficiency at an arbitrary 
range. 


INTRODUCTION 


The problem of microwave beam power transmission between two aperture antennas 
have been studied extensively (for example refs. 1-5) in open literature. The 
purpose of this paper is to determine the effective range for reliable and 
efficient beam power transfer between two aperture antennas. The transmitting 
tangential aperture field is simulated with a tapered aperture field distribu- 
tion. A numerical evaluation of Kirchoff and HelmOtz vector diffraction (ref. 
6) solution is used to calculate the electric and magnetic fields distribution 
at arbitrary observation distances. To obtain the transmission efficiency the 
power transmitted is compared to the power intercepted by receiving aperture 
antenna located at an arbitrary distance from the transmitting antenna (see 
Figure 1). A description and a copy of the program are included in appendixes 
A and B, respectively. 



POWER TRANSFER EFFICIENCY 


/ 


Let antennas A, and A 2 have generally dissimilar circular apertures of radii aj 
and a 2 , respectively, as shown in Figure 1. The antennas are separated by a 
length d along a common centerline (z-axis). This axis is perpendicular to the 
aperture planes of Aj and A 2 . The antennas are assumed to be matched to their 
feeding lines. Let assume that the transmitting antenna radiates an amount of 
power P t and the receiving antenna intercepts an amount of power P r . The 
transmission efficiency is given as: 

v ■ p r / P t (1) 

The transmitted power P t can be expressed as: 


P, = .5 // RE (E^x.y.O) x H,*(x,y,0)) ( z dS (2) 


and the received power P p can be expressed as: 


P r = .5 ff RE (E r (x,y,d) x H r *(x,y,d)) e z dS (3) 


where the asterisk indicates the complex conjugate values, e is the unit vector 
in the z direction, E t (x,y,d) is the tangential electric field component, and 
H t (x,y,d) is the tangential magnetic component of the field at an arbitrary 
location. The components H t (x,y,0), E r (x,y,d), and H (x,y,d) can be expressed 
in terms of E t (x,y,0) which uniquely dfetermines the field in the entire half- 
space z > 0 (ref. 7). The electric and magnetic fields at z * 0 are assumed to 
have the following forms respectively: 


E,(x,y,0) - E„(C + (1 - C)(l - (p / R) 2 ))" x (4) 
H,(x,y,0) = (e z x E,(x,y,0)) / Z 0 y (5) 

where the parameter C, p, R are defined as: 

C = 1(T et/20 , ET: edge taper in dB (6) 

p = x 2 + y 2 , radius (7) 

R , radius of aperture A 2 (8) 


The E r and H r fields due to the above aperture distribution are given by: 
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E r (x,y,d) = Jf E t (x,y,0) e' jkr ((l/r)(jK + 1/r) n. r) dS (9) 

H r (x,y,d) = (e z X E r (x,y,d)) / Z Q (10) 

Figure 2 depicts all required parameters to evaluate equations (9) and (10). 
Since the observation points are in the near field of the aperture, equations ( 9 ) 
and (10) are general, do not simplify further. In order to compute the fields 
numerically, we divide the aperture into a finite number of grid points and the 
total field at the n th observation point is the cartesian sum of the fields due 
to all grids on the transmitting aperture. 

A computer program based on equation (9) and (10) was developed to compute the 
transmitting electric field and power flow at an arbitrary distance. The program 
listing and a sample run is included in Appendix I and II respectively. 


DISCUSSION AND RESULTS 


The effects of range and aperture tapering were studied for two identical 
aperture antennas. The apertures and relevant parameters are described as 
fol lows: 


Transmitting aperture 
Receiving aperture 
Operating frequency 
Observation range 
Far field range (FF) 
Taper range 


A, = 0.7854 m 2 (Dia = 1 m) 

A ? = 0.7854 m 2 (Dia = 1 m) 
f = 3 GHz 

.1 FF < R < 1 FF (FF = 2D 2 /lam) 
FF = 20 m 

ODb , - lOdB , - 15dB , -20dB , -25dB , -30dB 


For the purpose of this report we will consider 50% efficiency to be the minimum 
in which beam power transmission will be useful. Figure 3 shows the resulting 
power transfer efficiency as a function of distance when the two apertures are 
of the same area. This results indicates that beam power is most effective at 
closer ranges. Uniform aperture illumination was used for the transmitting 
antenna. Figure 4 and show the results obtained for a cases when the transmit- 
ting antenna is two times larger than the receiving aperture. Figure 5 shows the 
results obtained for cases when the receiving antenna is two times larger than 
the transmitting antenna. Basically this result shows by increasing the 
transmitting antenna size does not necessary increase the transmission 
efficiency. In the contrary increasing the receiving antenna size will increase 
the amount of power that can be collect. Figure 6 shows the effect of tapering 
the transmitting aperture on the beam power efficiency at a given distance. 
Figure 5 indicates a greater efficiency number can be obtained by properly 
tapering the transmitting aperture. Figure 6 and 7 shows the previous cases 
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(figure 3 and 4) but with a -20 db taper level introduced in the transmitting 
aperture. In this case the useful range was increased by 20% . 

In general when we consider beam power transmission we have to carefully consider 
two parameters, the transmitting aperture size and the taper required to 
efficiently transmit power from point to point. It is also assumed in our 
discussion that the received antennas will be mechanically steered, therefore 
scanning will not be an issue. The results presented indicates that beam power 
should be considered as an alternative for point to point power distribution. 
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APPENDIX A 

DESCRIPTION OF PROGRAM 

A computer program was designed to calculate the power efficiency between two 
circular apertures at an arbitrary range. The method of analysis is vector-wave 
diffraction approach. The main features in this program (BEAM FORTRAN) are the 
calculation of power transfer efficiency, transmit power density and receive 
power density. The inputs are modified according to the case inside the program 
BEAM FORTRAN. The following is a description of all the inputs parameters to the 
the program BEAM. 


Input Parameters (BEAM FORTRAN) 


LAMBDA 

LF 

TAREA 
RAREA 
TOTAREA 
mts square 
DZ 
ET 


wavelenght in mts. 

spacing of computation grids 

transmit area in mts square 

receive area in mts square 

total area for calculating power radiated in 

fraction of far-field distance (2D**2/ lambda) 
edge taper of transmit aperture fields 


Example of Inputs and Results 

INPUTS: 


LAMBDA = . 1 

LF ' « .5 

TAREA * 0.7854 

RAREA = 0.7854 

TOTAREA - 1.5708 

DZ * .1 

ET = 0. 

RESULTS: 

FRACTION OF 2*D**2/LAMBDA = 0.100000024 
DISTANCE BETWEEN APERTURES = 0. 157079869E-02 KM 
EDGE TAPER = 0. 000000000 E+00 DB 

TRANSMITTING APERTURE AREA « 0.785399973 SQ. MTS 
RECEIVING APERTURE AREA = 0.785399973 SQ. MTS 
CAPTURE APERTURE AREA * 1.57079983 SQ. MTS 
TOTAL POWER RADIATED = 0.243498638E-03 W 
TOTAL POWER RECEIVED = 0.221619441E-03 W 
EFFICIENCY PERCENT = 91.0146332 % 
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APPENDIX B 


COMPUTER PROGRAM 

PROGRAM BEAM 

******* PROGRAM DESCRIPTION ***** 

* THIS PROGRAM CALCULATES THE BEAM POWER EFFICIENCY OF TWO APERTURES 

* AT AN ARBITRARY RANGE. THE INPUTS ARE THE WAVELENGTH (LAMBDA) IN 

* METERS, RECEIVING AND TRANSMITTING GRID SPACING AND AREA IN SQ MTS., 

* EDGE TAPER IN DBS AND THE SEPARATION BETWEEN THE APERTURES AS A 

* FRACTION OF 2D**2/ LAMBDA. 

******* VARIABLES ***** 

***** 


LF 

***** lambda 

***** RAREA 

***** tarea 
*****totarea 

***** Q2 
***** ET 
***** |< 
***** NO 
*****NXO NXO 
*****NXL NXL 


- GRID SPACING 

- WAVELENGTH 

- RECEIVING AREA 

- TRANSMITTING AREA 

- paptiirf area 

- FRACTION OF 2*D**2/ LAMBDA 

- EDGE TAPER 

- WAVE NUMBER 

- FREE SPACE IMPEDANCE 

- NUMBER OF POINTS IN RECEIVING GRID 

- NUMBER OF POINTS IN TRANSMITTING GRID 

- SQRT(-l) 

- EQUATION CONSTANT 

- FAR FIELD DISTANCE 

- EFFICIENCY 

- EFFICIENCY PERCENT 

- DELTAS FOR RECEIVING GRID 

- DELTAS FOR TRANSMITTING GRID 

- POINT IN RECEIVING GRID 

- POINT IN TRANSMITTING GRID 

- POWER TRANSMITTED 

- POWER RECEIVED 

- TRANSMITTED ELECTRIC FIELD 

- RECEIVED ELECTRIC FIELD 

- RECEIVED MAGNETIC FIELD 
******* VARIABLE DECLARATION ***** 

INTEGER ANO, BNO, CNL, DNL 
REAL LAMBDA, LF,M,KEM,K, NO 

COMPLEX ERX (300 , 300) , HRY (300 , 300) , ETX (300 , 300) , TREX , TRMY , J 
**************************** 

***** ***** input DATA ***** 


***** 

***** 


***** 

***** 


***** 


j 

KEM 
FF 
EFF 
PERC 
& DXY 
& DXY 


* DXY 

* DXY 
*****(XO,YO) 
*****(XL YL) 

***** ppi 

c**** 

***** 

* TREX 

* TRMY 


PR2 

ETX 

ERX 

HRY 


**************************** 

LAMBDA = .1 

LF * .5 

TAREA = 0.7854 

RAREA = 0.7854 

TOTAREA = 1.5708 

DZ = .1 

ET = 0. 

C**** *** CONSTANTS *** 

J = (0..1.) 

PI = 4 * ATAN(1.) 

NO = 120 * PI 

XAL = SQRT (TAREA) / 2 

XAO = SQRT (RAREA) / 2 

XAT = SQRT (TOTAREA) / 2 

K = 2. * PI / LAMBDA 

KEM = 1. / (4. * PI) 

DXY = LAMBDA * LF 

NXL = 2. * XAL / DXY + 1.1 

NXO = 2. * XAO / DXY + 1.1 


BEA00010 

BEA00020 

BEA00030 

BEA00040 

BEA00050 

BEA00060 

BEA00070 

BEA00080 

BEA00090 

BEA00100 

BEA00110 

BEA00120 

BEA00130 

BEA00140 

BEA00150 

BEA00160 

BEA00170 

BEA00180 

BEA00190 

BEA00200 

BEA00210 

BEA00220 

BEA00230 

BEA00240 

BEA00250 

BEA00260 

BEA00270 

BEA00280 

BEA00290 

BEA00300 

BEA00310 

BEA00320 

BEA00330 

BEA00340 

BEA00350 

BEA00360 

BEA00370 

BEA00380 

BEA00390 

BEA00400 

BEA00410 

BEA00420 

BEA00430 

BEA00440 

BEA00450 

BEA00460 

BEA00470 

BEA00480 

BEA00490 

BEA00500 

BEA00510 

BEA00520 

BEA00530 

BEA00540 

BEA00550 

BEA00560 

BEA00570 

BEA00580 

BEA00590 

BEA00600 
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NXT = 2. * XAT / DXY + 1.1 

M = 10**(-ET / 20) 

FF = (2. * (2. * XAL)**2. /LAMBDA) * DZ 

c**** 

c**** 

c**** 

C**** *** LOOP FOR TRANSMITTING GRIDS *** 

DO 50 CNL = 1, NXL 

XL = -XAL + DXY * (CNL - 1) 

DO 25 DNL = 1, NXL 

YL = -XAL + DXY * (DNL - 1) 

C**** *** TRANSMITTED ELECTRIC FIELD *** 

= M + (1. - M) * (1. - (XL**2 + YL**2)/ 
(XAL**2))**2 


GRIDS *** 


(ANO - 1) 

DO 100 BNO « 1, NXT 

YO = -XAT + DXY * (BNO - 1) 

Q**** 

C**** ***** SUBROUTINE FOR RECEIVED FIELDS CALCULATION ***** 

Q**** 

£**** 

CALL EANDH (M , K , J , NXL , XAL , DX Y , FF , NO , KEM , 

1 TREX , TRMY f ETX , XO , YO) 

ERX(ANO.BNO) = TREX 
HRY(ANO,BNO) = TRMY 
TREX =0 
TRMY = 0 

100 CONTINUE 

200 CONTINUE 


1 


ETX (CNL, DNL) 


25 

50 

0**** 

Qieidrk 

c**** 

£**** 

Q**** 


CONTINUE 

CONTINUE 


*** LOOP FOR RECEIVING 


DO 200 ANO = 1, NXT 
XO = -XAT + DXY * 


c**** 

0**** ***** SUBROUTINE FOR POWER AND EFFICIENCY CALCULATION ***** 

Q**** 

CALL FINISH (ERX,HRY,M, NXL, NXO.ET, ETX, KEM, LF,RAREA,TAREA, 

1 TOTAREA , XAL , XAO , DX Y , DZ , XO , YO , NO , FF , NXT , XAT) 

STOP 

END 


BEA00610 

BEA00620 

BEA00630 

BEA00640 

BEA00650 

BEA00660 

BEA00670 

BEA00680 

BEA00690 

BEA00700 

BEA00710 

BEA00720 

BEA00730 

BEA00740 

BEA00750 

BEA00760 

BEA00770 

BEA00780 

BEA00790 

BEA00800 

BEA00810 

BEA00820 

BEA00830 

BEA00840 

BEA00850 

BEA00860 

BEA00870 

BEA00880 

BEA00890 

BEA00900 

BEA00910 

BEA00920 

BEA00930 

BEA00940 

BEA00950 

BEA00960 

BEA00970 

BEA00980 

BEA00990 

BEA01000 

BEA01010 

BEA01020 

BEA01030 

BEA01040 

BEA01050 

BEA01060 

BEA01070 

BEA01080 

BEA01090 

BEA01100 

BEA01110 

BEA01120 

BEA01130 

BEA01140 


0**** ******************************************.************************g^Qj25Q 


C**** BEA01160 
0**** THIS SUBROUTINE CALCULATES THE XYZ COMPONENTS OF THE RECEIVED BEA01170 
0**** ELECTRIC AND MAGNETIC FIELDS FOR EVERY POINT IN THE RECEIVING BEA01180 
C**** APERTURE BY ADDING THE EFFECT OF EVERY POINT IN THE TRANSMITTING BEA01190 
0**** APERTURE. BEA01200 
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c**** 

0 **** 


0 **** 


Q'kick'k 

0 **** 


0**** 

0**** 

0**** 

0**** 

0 **** 

0**** 

0**** 

0***rt 

0**** 

0**** 

0**** 


0**** 

0**** 

0 **** 

0**** 

0**** 


0**** 

0**** 

0**** 

0**** 

0**** 

0**** 

0**** 

0**** 

0**** 

0**** 

0*rt*rt 

0**** 

0**** 

0**** 

100 

200 


BEA01210 

******************************************************************g£AQ1220 
SUBROUTINE EANDH (M, K , J , NXL , XAL , DXY , FF , NO, KEM , TREX , TRMY , ETX , XO , YO) BEA01230 


REAL M, KEM, K, NO BEA01240 

INTEGER CNL.DNL BEA01250 

COMPLEX J , PHASOR , EFOX , ESOX , ETX ( 300 , 300) , REX , RMY , TREX , TRMY BEAO 1 260 

*** LOOP FOR TRANSMITTING GRIDS *** BEA01270 

DO 200 CNL = 1, NXL BEA01280 

XL = -XAL + DXY * (CNL-1) BEA01290 

DO 100 DNL = 1, NXL BEA01300 

YL = -XAL + DXY * (DNL - 1) BEA01310 

BEA01320 

BEA01330 

*** DISTANCE *** BEAO 1340 

BEA01350 

BEA01360 

BEA01370 

R = SQRT((X0-XL)**2 + (Y0-YL)**2 + FF**2) BEA01380 

BEAO 13 90 

*** COS X *** BEA01400 

BEA01410 

BEA01420 

COSX = FF / R BEA01430 

BEAO 1440 
BEA01450 
BEA01460 

*** FIRST ORDER RECEIVED ELECTRIC FIELD *** BEA01470 


BEA01480 
BEA01490 
BEA01500 
BEA01510 
BEA01520 
BEA01530 
BEA01540 
BEA01550 
BEA01560 
BEA01570 
BEA01580 
BEA01590 
BEA01600 
BEA01610 
BEA01620 
BEAO 1630 
BEA01640 
BEA01650 
BEA01660 
BEA01670 
BEA01680 
BEA01690 
BEA01700 
BEA01710 
BEA01720 
BEA01730 
BEAO 1740 
BEAO 17 50 

TREX = TREX + REX BEA01760 

TRMY * TRMY + RMY BEA01770 

CONTINUE BEA01780 

BEA01790 

CONTINUE BEA01800 


PHASOR = (COS(K*R) - J*SIN(K*R>) / R 

EFOX = J * K * ETX (CNL, DNL) * COSX * PHASOR 


*** SECOND ORDER RECEIVED ELECTRIC FIELD *** 


PHASOR * (COS(K*R) - J*SIN(K*R)) / R**2 
ESOX = ETX (CNL, DNL) * COSX * PHASOR 


*** COMBINED RECEIVED ELECTRIC ELECTRIC FIELDS *** 


REX = KEM * (EFOX + ESOX) * DXY * DXY 
*** RECEIVED MAGNETIC FIELD *** 

RMY = REX / NO 


*** TOTAL FIELDS *** 
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RETURN BEA01810 

END BEA01820 

q**** ******************************************************************g£^Q1830 


C**** BEA01840 
C**** THIS SUBROUTINE CALCULATES THE TOTAL POWER RECEIVED, TOTAL BEA01850 
C**** • BEA01860 
C**** powER TRANSMITTED AND EFFICIENCY PERCENT. BEA01870 
C**** BEA01880 


0**** ******************************************************************B£AQ}ggQ 


Q**** 

c**** 

SUBROUTINE FINISH (ERX, HRY, M,NXL,NXO, ET,ETX,KEM,LF,RAREA,TAREA, 
1 TOTAREA,XAL,XAO,DXY,DZ,XO,YO,NO,FF,NXT,XAT) 

REAL LF,M, KEM,NO 
INTEGER ANO.BNO 

COMPLEX ERX (300 , 300) , HRY (300 , 300) , ETX (300 , 300) 

0**** 

c **** ***** T0T AL POWER RADIATED ***** 

0**** 

PR1 = 0 

DSZO = DXY * DXY 
DO 200 ANO = 1, NXT 

DO 100 BNO = 1, NXT 
DSZO = DXY * DXY 

PR1 = .5 * REAL (ERX (ANO, BNO) * CONJG (HRY (ANO , BNO) ) ) 

1 * DSZO + PR1 

100 CONTINUE 

200 CONTINUE 
0**** 

0**** 

0 **** 

0**** ***** TOTAL POWER RECEIVED ***** 

0**** 

PR2 = 0 

N1 = (-XAO + XAT) / DXY + 1.1 
N2 = (XAO + XAT) / DXY + 1.1 
DO 800 ANO = Nl, N2 

DO 900 BNO = Nl, N2 

PR2 = .5 * REAL (ERX (ANO, BNO) * CONJG (HRY (ANO, BNO))) 

1 * DSZO + PR2 

900 CONTINUE 

800 CONTINUE 

0**** ***** EFICIENCY PERCENT CALCULATION ***** 

0 **** 

EFF = PR2 / PR1 
PERC = EFF * 100 

0**** 

0**** *** OUTPUT DATA *** 

WRITE(9,*) 'FRACTION OF 2*D**2/LAMBDA = ',DZ 
WRITE(9,*) 'DISTANCE BETWEEN APERTURES = ' ,FF/1000, ' KM' 
WRITE(9,*) 'EDGE TAPER = * , ET , ' DB ' 

WRITE(9,*) 'TRANSMITTING APERTURE AREA = ' ,TAREA, ’ SQ. MTS' 
WRITE(9,*) 'RECEIVING APERTURE AREA = ' ,RAREA, 'SQ. MTS' 
WRITE(9,*) 'CAPTURE APERTURE AREA = ' .TOTAREA, 'SQ. MTS' 
WRITE(9,*) 'TOTAL POWER RADIATED = ' , PR1 , 'W* 

WRITE(9,*) 'TOTAL POWER RECEIVED = \PR2,'W' 

WRITE(9,*) 'EFFICIENCY PERCENT = ' ,PERC f '%' 

RETURN 


BEA01900 

BEA01910 

BEA01920 

BEA01930 

BEA01940 

BEA01950 

BEA01960 

BEA01970 

BEA01980 

BEA01990 

BEA02000 

BEA02010 

BEA02020 

BEA02030 

BEA02040 

BEA02050 

BEA02060 

BEA02070 

BEA02080 

BEA02090 

BEA02100 

BEA02110 

BEA02120 

BEA02130 

BEA02140 

BEA02150 

BEA02160 

BEA02170 

BEA02180 

BEA02190 

BEA02200 

BEA02210 

BEA02220 

BEA02230 

BEA02240 

BEA02250 

BEA02260 

BEA02270 

BEA02280 

BEA02290 

BEA02300 

BEA02310 

BEA02320 

BEA02330 

BEA02340 

BEA02350 

BEA02360 

BEA02370 

BEA02380 

BEA02390 

BEA02400 
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END 


BEA02410 
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/*EXEC TO RUN BPE*/ 

SETUP FTN 
"FTNLIB" 

"FI 09 DISK CAACOSTA BPE al" 
"LOAD beam (start" 
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Figure 1. Geometry for calculating 
the transmission efficiency 



Figure 2. Parameters description for obtaining 
beam power efficiency 
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EFFICIENCY PERCENT EFFICIENCY PERCENT 



Figure 3. Efficiency curves resulting from equal areas 
with two different tapers 



Figure 4. Efficiency curves resulting from different 
transmitting - and receiving - aDertures size . 
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EFFICIENCY PERCENT EFFICIENCY PERCENT 



DISTANCE (FRACTION OF 2D**2/LAMBDA) 

Figure 6. Efficiency curves resulting from different 
transmitting and receiving apertures size 
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Y PERCENT EFFICIENCY PERCENT 



0 1 Z 3 4 5 

RECEIVING APERTURE DIAMETER (METERS) 


Figure 7. Efficiency curves for different receiving 
aperture size at far-field range 



Figure 8. Efficiency curve for different transmitting 
aperture tapers at near-field range. 
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